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Abstract We discuss the concept of the Berry phase in a dissipative system. We 
show that one can identify a Berry phase in a weakly-dissipative system 
and find the respective correction to this quantity, induced by the envi- 
ronment. This correction is expressed in terms of the symmetrized noise 
power and is therefore insensitive to the nature of the noise representing 
the environment, namely whether it is classical or quantum mechani- 
cal. It is only the spectrum of the noise which counts. We analyze a 
model of a spin-half (qubit) anisotropically coupled to its environment 
and explicitly show the coincidence between the effect of a quantum 
environment and a classical one. 
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Introduction 

Three papers published independently in 1932 by Zener, Landau and 
Stueckelberg (Landau, 1932; Zener, 1932; Stueckelberg, 1932) have in- 
troduced the phenomenon known today as Landau-Zener tunneling. The 
idea is to consider a 2-level system, where the energy of each level varies 
linearly with a classical variable (which, in turn, is varied linearly in 
time). As function of time, t, the energy levels should intersect but for 
the inter-level coupling A which gives rise to an "avoided crossing" in the 
spectrum, cf. Fig.l. Using the spin notation, one can write the Hamil- 
tonian as 7^ = atS^ + ASx- Here S = a/2, and cr^, ax are Pauli spin- 1/2 
operators; a is the rate of change of the energy of the pseudo-spin at 
asymptotic times. The avoided crossing gap is A. The probability of 
transition from, say, the lower level at time — oo , to the upper level at 
time -l-oo is given by Plz = exp[— (7r/2)A^/a!]. 

Besides being ubiquitous in physics and chemistry, the Landau-Zener 
framework appears to suggest a natural definition for the notion of 
adiabaticity. The adiabatic limit is approached when i^z << 1, i-e., 
a « A^. The latter inequality involves a comparison of the rate of 
change (of the time dependent term in the Hamiltonian) with the gap 
in the spectrum, A. This notion of the adiabatic limit has become 
widespread. A closer look suggests that, in general, adiabaticity cannot 
be associated with comparing the rate of change to the gap. Indeed, 
on one hand any finite, discrete-spectrum system is coupled, however 
weakly, to the rest of the universe. Hence the emerging spectrum is, at 
least in principle, always continuous and gapless. The naive view would 
then imply that the adiabatic limit cannot be approached. This, on the 
other hand cannot be correct: if we consider a finite system with a dis- 
crete spectrum, for which adiabaticity is well defined, it is inconceivable 
that an infinitesimal coupling to the continuum (rendering the overall 
spectrum continuous) will change its physics in a dramatic way. The 
resolution of this problem is provided by the observation that the cri- 
terion for adiabaticity involves not only spectral properties but also the 
matrix elements of the system-environment coupling. 

To gain some insight into this problem we focus here on the analysis 
of the Berry phase (Berry, 1984) in a weakly dissipative system. It is 
particularly timely to address this issue now given the recent experi- 
mental activities in realization of controlled quantum two-level systems 
(qubits), and in particular, the interest in observing a Berry phase (BP) 
(see, e.g., (Falci et al., 2000)). For instance, the superconducting qubits 
have a coupling to their environment, which is weak but not negligi- 
ble (Nakamura et al., 1999; Vion et al., 2002; Chiorescu et al., 2002), 
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and thus it is important to find both the conditions under which the 
Berry phase can be observed and the nature of that Berry phase. 

In this paper we appeal to a simple analysis of the problem. We 
first, in Section 2, consider a quantum-mechanical framework, where a 
perturbative approach is taken. When the environment is replaced by a 
single oscillator, a second-order perturbation analysis is straightforward 
and produces a result which allows for a simple interpretation. We then 
generalize the calculation for a host of environmental modes. In Section 3 
we consider a toy model where the environment is replaced by a classical 
stochastic force. The quantities of interest, the Lamb shift and the Berry 
phase, are then calculated, and simple heuristic arguments are given to 
interpret the results. To complete the analogy with the analysis of the 
previous section, here the "single-oscillator environment" is replaced by 
a simple periodic classical force (of random amplitude). In Section 4 we 
summarize the relation between the quantum mechanical approach and 
the classical model in more general terms. 

1. The system: spin + environment 

We begin in the conventional way by writing the Hamiltonian for the 
"universe" (system -|- environment) as 



The system is defined as the set of those quantum degrees of freedom 
that one is interested to control and measure; the environment consists of 
all the rest, namely those degrees-of-freedom we can neither control nor 
measure. The coupling between the system and environment is Vcoupiing- 
The properties of the environment are controlled by macroscopic param- 
eters, such as temperature. Our treatment below applies to a reservoir 
at either zero or a finite temperature. 

For our purposes it is sufficient to represent the environment by a sin- 
gle operator X which couples to a spin. The Hamiltonian then becomes 



Hereafter we put ng = 1. Below we express our results in terms of the 

statistical properties (correlators) of the environment's noise, X{t). De- 
pending on the physical situation at hand, one can choose to model the 
environment via a bath of harmonic oscillators (Fcynman and Vernon, 
1963; Caldeira and Leggett, 1983). In this case the generalized coordi- 
nate of the reservoir is defined as X = Y1, ^i^i, where {xi} are the coor- 
dinate operators of the oscillators and {Aj} are the respective couplings. 
Eq. 2 is then referred to as the spin-boson Hamiltonian (Leggett et al.. 




(1) 



n = -\ligB-&- \Xu^ + ii, 
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1987). Another example of a reservoir could be a spin bath (Prokof'ev 
and Stamp, 2000) ^. However, in our analysis below we do not specify 
the type of the environment. Wc will only assume that the reservoir 
gives rise to markovian evolution on the time scales of interest. More 
specifically, the evolution is markovian at time scales longer than a cer- 
tain characteristic time Tc, determined by the environment ^. We assume 
that Tc is shorter than the dissipativc time scales introduced by the envi- 
ronment, such as the dephasing or relaxation times and the inverse Lamb 
shift (the scale of the shortest of which we denote as Tdiss, Tc <C Tdiss)- 
We further assume that Tc <^ tp, the characteristic variation time of 
the field B(t). Moreover, under these conditions we may consider only 
lowest-order (in the system-environment coupling) contributions to the 
quantities of interest: energy shifts, BP and relaxation rates. Indeed, 
if one divides the evolution time interval into short domains (<C tp), 
longer than Tc but shorter than Tjiss, fluctuations at different domains 
are uncorrelated and their effect can be analyzed separately. At the same 
time, for each domain (^ Tdiss) the effect of the noise is weak. Thus, 
to the leading order corrections to the dynamics may be described as 
corrections to the rates (energies) of the spin dynamics, which may be 
estimated perturbatively. We also consider an underdamped spin, with 
the dissipativc times longer than the period of the coherent dynamics, 
Tdiss ^ 1/B. This implies that the time windows alluded to above 
consist of numerous oscillations, in other words they are S> 1/-B. 

We have chosen an anisotropic spin-environment coupling, oc a^. This 
is a realistic model, e.g., for many designs of solid-state qubits, where 
the different components of the "spin" are influenced by entirely different 
environmental degrees of freedom (Nakamura et al., 1999; Vion et al., 
2002; Chiorescu et al., 2002). While our analysis can be generalized to 
account for multiple-directional fluctuating fields (Whitney et al., 2004), 
here we focus on unidirectional fluctuations (along the z axis). 



^For any reservoir in equilibrium the fluctuation-dissipation theorem provides the rela- 
tion between the symmetrized and antisymmetrizcd correlators of the noise: Sxi'^) = 
Ax{i^) coth(w/2T). Yet, the temperature dependence of Sx and Ax may vary depending 
on the type of the environment. For an oscillator bath, Ax (also called the spectral den- 
sity Jx{<^)) is temperature-independent, so that Sx{^) = Jx{^) coth(to/2T). On the other 
hand, for a spin bath Sx is temperature-independent and is related to the spins' density of 
states, while Ax{i^) = Sjf (ti^) tanh(a;/2T). 

^This time may be given by the correlation time of the fluctuations, but in general is a more 
subtle characteristic of the spectrum related to its roughness near qubit's frequencies. Note 
further that for singular spectra Tc may be ill defined and the perturbativc analysis may fail. 
See, e.g., (Bloch, 1957; Redfield, 1957; Slichter, 1978; Makhlin et al., 2003; Wilhelm et al., 
2004; Whitney and Gefen, 2004). 
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Another remark to be made concerns the possibility to observe a 
(weak) dissipative correction to Berry phase in spite of the dephasing 
and relaxation phenomena. While the respective time scales (Ti, T2 and 
the inverse of the correction to the Berry phase) scale similarly with the 
strength of fluctuations (inversely proportionally to the noise power), 
they are dominated by different frequency domains. Indeed, the dephas- 
ing and relaxation arc known to be dominated by resonant fluctuations 
with frequencies close to B (for the relaxation and the corresponding 
contribution to dephasing) and (for the pure dephasing), cf. Eq. (15) 
below. In contrast, as we shall see below, the Lamb shift and the correc- 
tion to the Berry phase accumulate contribution from the entire range 
of frequencies. Thus, one may think of (engineering) a system with an 
environment whose fluctuations at ~ i? and 1/ ~ are suppressed. In 
this case, one can easily observe an observable correction to the Berry 
phase at times when the dephasing and relaxation are still negligible. 

2. Quantum- mechanical analysis 

In this section wc consider a two-level system coupled to an environ- 
ment which we treat as a quantum-mechanical system. We begin with a 
discussion of the Lamb shift and then show, in Subsection 2.3, how the 
results for the Lamb shift may be used to find the environment-induced 
correction to the Berry phase and the relaxation times. 

2.1 Lamb shift as level repulsion 

Consider first, for illustration, a simple system of the spin coupled to 
a single oscillator, with the Hamiltonian 

H = —^Baz — ^cax{a^ + a) + ujQa^a , (3) 

where c is the coupling constant. Let \n) denote the n-th level of the 
oscillator; the second-order corrections to the energies of the states |t, 0) 
and 11,0) are 

^(2)_ |(T,o|v|i,i)|^ _ 1 

T Luo + B 4:LUo + B' ^ ' 

and 

^p)^JU.o|v|T,i)l^_i^^ 

^ UQ-B AujQ-B ^ ' 

where V = {c/2) Gx {o) + a) is the perturbation. This results in the 
following correction to the level spacing E\ ~ Ey. 



„2 
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This correction (the Lamb shift) has different signs for fast (luq > B) 
and slow {ooq < B) oscillators. As one can see from Eqs. (4), (5), this 
result can be understood in terms of the level repulsion (Wilhelm et al., 
2004): the perturbation couples the level |T, 0) to |i, 1) and |i, 0) to |T, 1). 
The levels of the latter pair are closer, and the coupling has a stronger 
effect on their energies. They repel each other due to the coupling, thus 
reducing the distance between |T, 0) and |J,, 0) for cjq > -B and increasing 
it for ujQ < B. 

2.2 Second-order perturbative analysis 

In this section we find the Lamb shift using the lowest-, second-order 
perturbative analysis. In the Hamiltonian (2) we treat the coupling term 
V = —\Xaz as a perturbation: H = Ho + V. The eigenstates of Hq 
are \a,i), where a =13/13 denotes the eigenstates of the spin without 
dissipation, with the spin direction parallel or antiparallel to the filed B, 
and i denotes eigenstates of the environment. The perturbation theory 
gives for the corrections to their eigenenergies: 

p(2) ^ _ V- \{a,i\V\f^,j)f (r,. 

h Ef + i^f - E^^^ - Ef^ - iO • 

For V = -^^Xgz we notice that (T_b| o-^ |Tb)^ = (iB|o"z|iB)^ = cos^ 6' 
and {'\b\ <^z IIb)^ = (isl '^z iTs)^ = sin^ 6*, and find for the environmcnt- 
averaged quantities Ea = pi E^ ^ (see the discussion of these quan- 
tities at the end of this subsection): 

(2) _ cos^6> ^ Pi I {i\X\j) 1^ _ sin^e ^ pi \ {i\ X \j) |^ 

' " 4 2_.^(o)_^(o)_.o 4 2.^^^(o)_^(o)_.o- 

(8) 

The correction to E]^ is obtained by substituting B — > —B into the above 
equation. Now using the identity 

oo 

^^g-i(B-iO)t^ (9) 



£;-iO 
we rewrite Eq. (8) as 

Ef^ = -^-j^dt {X{t)X{Q)) (cos^ e + sin^ ^e-^^^*) 6"°* , (10) 

where we have used the relation 

{X{t)Xm =^p, {i\ X \j) {j\ X |z) e-(^^-^')* . (11) 
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In terms of the the Fourier transform {X'^) = J dt {X[t)X{0)) e*'^* we 
obtain 

= os=« f^Ji (f).^. (12) 

T 4 J 27r - iO 4 J 27r + S - iO ^ ^ 

For the Lamb shift E^^^^ = ^e{E^^^ — -Ej^-*) this gives a principal value 
integral 

^La^b - 2 Sin --^ - 5sm Vj^ - -^rzr^, , (13) 
where 

Sx{v)^\{{Xl) + {Xl,)) = \j dt([X(i),X(0)]+)e-*. (14) 

Thus the Lamb shift is expressed in terms of the symmetrized correlator 
Sx and is insensitive to the antisymmetric part of the noise spectrum. 

As one can see from Eq. (13), in agreement with the discussion in 
the previous section, the high-frequency noise [u > B) reduces the en- 
ergy gap between the spin states (Leggett et al., 1987), while the low 
frequency modes {v < B) increase the energy gap. 

Similarly, from Eq. (12) one can evaluate the dephasing time: 

1 = - am(£f > + Ef^) = "^S.iu = 0) + '^S,{u = B) . (15) 

This expression correctly reproduces the contribution of the transverse 
fluctuations (oc sin^ 6) to the dephasing rate, but underestimates the 
longitudinal contribution (oc cos^ 9) by a factor of two (cf. Rcf. (Bloch, 
1957; Redfield, 1957; Weiss, 1999)). One can show that an accurate 
evaluation of this contribution, as well as the analysis of the relaxation, 
requires taking into account corrections to the eigenstates, and not only 
to the eigenenergies (7). More precisely, our calculation of the correc- 
tions to the eigenenergies in this subsection corresponds to evaluation 
only of the four left diagrams in Fig. 7 of Ref. (Makhlin et al., 2003); 
the term iO in the denominators allows one to find also the outgoing 
transition rates from the eigenstates (and the respective contribution, 
oc sin^ 9, to dephasing) but only the part of the 'pure-dephasing' rate, 
oc I cos^ 9. Analysis of the two remaining diagrams in Fig. 7 and those 
in Fig. 6 allows one to find also the pure dephasing rate (as well as 
the incoming transition rates, the latter though do not require an extra 
evaluation due to probability conservation). 
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2.3 From Lamb shift to Berry phase 

So far we have analyzed the environment-induced correction to the 
level splitting (the Lamb shift). Using the results above one can eval- 
uate also the environment-induced correction to the Berry phase for 
a slow cyclic variation of the magnetic field B (Whitney and Gefen, 
2001; Whitney and Gefen, 2003; Whitney et al., 2004; Whitney and 
Gefen, 2004). 

Indeed, consider the simplest case of conic variations of the field 
around the z-axis (to which the environment is coupled), as shown in 
Fig. 1: the field varies at a constant rate, with the low angular velocity 
ujb, and traverses the circle after the period tp = 27r/c<jB- The analysis 
of the spin dynamics is considerably simplified by going to the frame, 
rotating with the angular velocity ujbz, where z is the unit vector along 
the 2-axis. In this frame the spin is subject to the fluctuating field Xz 
and the field B + lo-qz, which is stationary. Thus, in this frame one can 
use the results of the analysis above to obtain the Lamb shift, if one sub- 
stitutes B by B + uj-qz. In other words, the correction to the Lamb shift 
associated with the variation of the field B in time, is given by taking the 
derivative oj-qQez of the Lamb shift (13) and multiplying by the period of 
variation, tp. After a full period the basis of the rotating frame makes a 
complete circle and returns to its initial position, i.e. coincides with the 
laboratory frame's basis. Hence the phases accumulated in the rotating 
and laboratory frames coincide, and it is sufficient to evaluate it in the 
rotating frame. Thus, one finds the environment-induced correction to 
the Berry phase to be 

5$BP = 27r^^^|^. (16) 



Taking the derivative of Eq. (13), we find: 

5(^)$BP = cos^sm^ e vj d. ^f^lgr.'f^ • (17) 

(Notice the convention: this expression gives the correction to the rela- 
tive Berry phase between the spin-up and spin-down states, rather than 
to the phases of each of these states.) As for the Lamb shift, the contri- 
butions of the high- and low-frequency fluctuations are of opposite signs. 
For the Berry phase the contribution changes sign &t v = 3B/2. 

In passing we note that this analysis can be generalized to an arbitrary 
(but adiabatic) path B(t), this enables one to see that the correction 
to the Berry phase is geometric, but that its geometric nature is very 
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Figure 1. 

different from the Berry phase of an isolated spin-half (Whitney et al., 
2004). 

In Section 3 we shall find exactly the same expression for the Lamb 
shift and therefore for the Berry phase in the case of classical environ- 
ment. 



2.4 High-frequency noise: renormalization of 
the transverse B-field 

Consider now the influence of the high-frequency fluctuations in the 
environment only [v ^ B). Since the frequencies of the fluctuations 
are much higher than the typical spin-dynamics frequencies, one may 
eliminate these high-frequency fluctuations using the adiabatic (Born- 
Oppenheimer) approximation, as described, e.g., by Leggett et al. (Leggett 
et al, 1987). 

Indeed, consider the spin-boson model, with the Hamiltonian 



H = -\{B + Xz)(j + n^ 



(18) 



where X = X^jCj(a| -|- a^) and Hem = J2i^i<^lc'i- Let us ignore the 
low- frequency oscillators and focus on those at high frequencies v ^ B. 
These fast oscillators adjust almost instantaneously to the slowly varying 
spin state. For the last two terms of the Hamiltonian (18) two lowest- 



energy states are 



T = IT)n^ gl) and i 



ii)n^ 



Here 



9i 



denotes the ground state of the ith oscillator corresponding to the spin 



state It), i.e. the ground state of uj-a-a - +c-{al + a^), and g^j is defined 
similarly; further eigenstates of the last two terms are separated by a 
gap ~ v. 
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Consider now the matrix elements of the first term —^Ba in this 
two-state low-energy subspace; one finds that its transverse component 
is suppressed by the factor 

ri(«ik.') =n»i>(-?/2-?) -xp(-/;|^) . (i9) 

where Jx(i^) = ti" Z^j c|(^(i^ — Ui) is the spectral density of the oscillator 
bath. At a finite temperature T each high-frequency oscillator remains 
in its thermal equilibrium state (subject to the spin state), rather than 
the ground state, and on the rhs of Eq. (19) the spectral density Jx{^) 
is replaced by the thermal noise power Sx{i^) = Jx{t^) coth(z//2fcBT')- 

Thus the role of the high-frequency oscillators is to suppress the trans- 
verse field component (in other words, the transverse ^f-factor). If we are 
interested only in the contribution to the level spacing (the Lamb shift), 
one should consider only the longitudinal (|| -B) part of the renormaliza- 
tion, i.e. multiply the result by sinO, to obtain Eq. (13). 

2.5 Effective-action analysis 

One can study the spin dynamics integrating out the environment and 
using the effective action for the spin. We derive the effective action us- 
ing the Feynman-Vernon-Kcldysh technique. For the interaction —Xsz 
with the z-component of the spin, the effective action (the influence 
functional) reads 

i$infl = -J/ dt [ dt's,{t)-s,{t')[iGx{t,t')], (20) 
2 JCk JCk 

where we assumed the Gaussian statistics of X, and defined the Green 
function Gx as iGxit^t') = {Tcj(X(t)X{t')). The time ordering here 
refers to the Keldysh time contour Ck, and in Eq. (20) we integrate 
over Ck', accordingly each of the time dependent variables assumes a 

'Koldysh index' u, d indicating the upper /lower branch of this contour. 
After the Keldysh rotation one obtains the influence functional in 

terms of the classical and quantum components, s§ = (s" + sf)/2 and 

o9 = _ od. 
— ^z ^z- 

^infl = - / dtdt' [sl{t)G^{t - t')sl(t') + \ sl{t)Gl{t - t')sl{t') 

(21) 

in terms of the retarded and Keldysh Green functions, G\ = —i6{t — 
t'){[X{t),Xit')U) and G| = -i([X(t), X(i')]+) = -2iSx{t-t'). 
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For classical noise X the commutator in the definition of vanishes, 
and one finds 

^h^r = ^ / / di' slit) Sx{t - t') slit') . (22) 

The results (13), (17) for the Lamb shift and the Berry phase involve 
only Sx and not the antisymmetrized correlator. Hence for the analysis 
of these quantities it should be sufficient to use the functional (22). 
Alternatively, one may consider a problem with a classical random field 
X{t) to reproduce these results. In the next section we perform the 
corresponding analysis. 

3. The Classical Model 

In this section we analyze the dynamics of a spin subject to a classical 
random field and derive the equation of motion for the spin dynamics 
(the spin-evolution operator), averaged over the fluctuations. Following 
the discussion of the case with quantum fluctuations, we first analyze the 
dynamics in a stationary field B and a random field; exactly as in the 
quantum case one can reduce the analysis of the dissipative corrections 
to the Berry phase accumulated over a conic loop to the problem with 
a stationary field by going over to a rotating frame. 

As we have demonstrated above, in the quantum problem the results 
for the corrections to the phase and dephasing, associated with the con- 
trolled dynamics of the magnetic field, involve only the symmetric part 
of the noise correlator, one expects that the results for these quantities 
in the classical problem, expressed in terms of the noise power, would 
coincide with the quantum results. Indeed, we find this relation below. 

Specifically, we analyze the following problem: a spin S is coupled 
to a controlled magnetic field B (stationary for now, but to be varied 
slowly in a Berry-phase experiment) and a randomly fluctuating field 
X{t), which we treat as a random variable with the correlation function 
given by Sx{t). Its dynamics is governed by the Larmor equation: 

S = [B + X{t)]x S . (23) 

This equation can be used to describe the dynamics of either a classical 
spin or the average spin value (i.e. the density matrix) of a spin-1/2. 

As we discussed in the Introduction, we assume that the noise is weak 
and short-correlated, i.e., that considerable dissipative contributions to 
the spin dynamics arise on time scales much longer than the typical 
correlation time of the noise. Below we discuss the influence of the 
low- and high-frequency fluctuations on the (classical) spin dynamics 
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and recover the results of the quantum analysis above. Further, using 
the result for the low-frequency contribution we obtain the correction to 
the Berry phase from the environmental fluctuations at all frequencies. 

3.1 Low- frequency noise: Lamb shift 

Consider first the effect of a slowly fluctuating random field X = Xz. 
Similar to the quantum-mechanical analysis in Section 2 we begin with 
the case of harmonic fluctuations (of random amplitude) and purely 
transverse noise {B = Bx, i.e. = 7r/2). Consider fluctuations X = 
Cjy cos{vt) at a low frequency v <^ B, during a time interval St. To 
evaluate the evolution operator, we analyze the dynamics in a reference 
frame (^,?},C) fluctuating together with the field (with the ^-axis along 
B + X{t) and the ry-axis, for instance, orthogonal to B and -X"). Since 
the fluctuating angular velocity of this frame's rotation is negligible, ~ 
Cyv/B <C Cy, the effective magnetic field in this frame \B + X{t)\C points 
along the (^-axis. Thus the dynamics reduces to rotation about this axis 
by the angle (t>{t) = //+'^* dr|S+X(r)| !^ //+'^* dr(S+X2(r)/2B), where 
B = \B\. Averaging the transverse spin component Sx + iSy oc e*''^^*^ 
one finds a lowest-order contribution to the phase factor, 5t{X'^)/2B, 
i.e. a Lamb shift c^/AB (we assumed 6t much longer than the period of 
oscillations, l/u). 

In principle, the evolution in the laboratory frame differs from that in 
the rotation frame. Transformation to/from the rotation frame at the 
beginning and the end of the time interval introduces corrections to the 
evolution operator or order c^/B. This is however a negligible boundary 
contribution. Indeed, for a sufficiently long time interval 5t ^ ^/cy the 
phase shift due to the Lamb shift, of order c^5t/B, is much larger (but 
still small, as long as 5t <^ B/c^). 

Similar results hold for more general low-frequency fluctuations, non- 
harmonic and with arbitrary direction 0. Indeed, in the same rotating 
frame the dynamics reduces to rotation about the i^-axis by the angle 
cp{t) = //+^* dT\B + X{t)\ ^ //+'^* dT{B + X\\{t) + Xi(r)/2S), where 
X\\ = X cosO, X± = X sinO are the longitudinal and transverse compo- 
nents of X (relative to B). Averaging the transverse spin component 
oc e^'^(*) one flnds, apart from dephasing, a lowest-order contribution to 
the phase factor, St{X'j)/2B, and hence the Lamb shift 

^^.,„^,/^lMd, (24) 

where 9 is the angle between B and the direction z of the noise. This 
result coincides with the low-frequency contribution in Eq. (13). 
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3.2 From low frequencies to all frequencies 

The expression (24) and the symmetry of the problem suggests a way 
to find the contribution of ah, not only slow, modes in the environment to 
the Lamb shift (and later to the Berry phase). Indeed, we discuss weak 
short-correlated noise, i.e. such that its contribution to the dynamics 
on time scales or order Tc is small. Contributions from different time 
intervals ~ Tc are uncorrelated and add up independently. Hence in 
the evaluation of the (real and imaginary) contribution of such a short 
interval to the evolution frequencies (the Lamb shift, the dephasing and 
relaxation rates) it is enough to consider the lowest nonvanishing, i.e. 
second order. 

The symmetry of the problem can be used to analyze the structure of 
such a second-order contribution. The spin-rotational symmetry (about 
the S-field's direction) and the time-translational symmetry imply that 
(i) the longitudinal and transverse fluctuations, X± and Xy, do not 
interfere and may be considered separately; (ii) it is convenient to expand 
the transverse fluctuating field in circularly polarized harmonic modes, 
and the latter contribute independently. 

The longitudinal noise gives rise to the pure dephasing (and only the 
low frequencies ^ I/T2 contribute), without affecting the level splitting. 
As for the transverse noise, for a single circularly polarized mode at fre- 
quency ly it is convenient to analyze its contribution in the spin frame, 
rotating at frequency v around the field B. In this frame the Larmor 
field is B — in the direction of B, and the fluctuating circularly polar- 
ized mode is slow. Applying to this mode Eq. (24), going back to the 
laboratory frame and adding up contributions of all modes we arrive at 
the expression for the correction to the Larmor frequency: 

SE = sin' evf ^4^. (25) 
J Air B — u 

It is thus this result which needs to be compared with the quantum 
correction (Lamb shift) of the previous section. Symmetrization of the 
integral in Eq. (25) w.r.t. to v brings it to the form of Eq. (13). Notice 
that regularization of this expression via the introduction of +iO in the 
denominator allows us also to recover the imaginary part of the Lamb 
shift, i.e. a contribution to the dephasing rate. 

3.3 High frequencies 

Although Eq. (25) describes the contribution of all frequencies, it 
is interesting to discuss specifically the limit of high frequencies. In 
this subsection we provide an argument which parallels the result of 
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subsection 2.4: the high-frequency fluctuations {u S> B) suppress the 

transverse (_L z) component of the S-field. 

Indeed, to solve for the dynamics in the presence of high-frequency 
fluctuations in a fixed direction, X{t)z, and the static field B, let us 
analyze the dynamics in the frame that rotates about the z-axis with 
angular velocity X{t), i.e. differs from the lab frame by a rotation by 
the fluctuating angle ^{t) = /QX(r)dr. The rotation of this frame is 
chosen to exactly compensate for the field X{t)z, and the Larmor field 
B(t) in this frame is just the S-field, but now fluctuating due to the 
frame's rotation as shown in Fig. 2. The spin dynamics is governed by 
the Larmor equation S = B{t) x S, and the value of the spin changes 
considerably only on time scales of order 1/B, during which many fluc- 
tuations occur. Looking at the dynamics on intermediate time scales, 
between l/v and 1/B, one finds that the spin dynamics is governed by 
the value of the S-field averaged over fast fluctuations. The averaging 
affects only the horizontal (orthogonal to z) component of the S-field. 
The direct evaluation shows that the horizontal component is suppressed 

oo 

exactly by the factor exp[- / {di'/2Tr)Sxii^)/'^'^] (cf. Eq. (19)). For in- 



stance, for a single mode at frequency ly we have X{t) = 2XjyCos(z/t) 
and <&(t) = 2Xys'm(i't)/i'; then the transverse component of the field is 
suppressed by the factor 1 - ($^)/2, and ($^)/2 = {XD/iP'. This evalu- 
ation of the dynamics in the rotating frame relies on the small parameter 
B/v. 

The spin-evolution operator (before averaging) Oiab(^)*') in the lab- 
oratory frame is related to that in the rotating frame, Oiab(i,i') = 
dz{-^{t))d^ot{t, t')dz{^{t')), via the transformation ©^(^(t)) from the 
lab frame to the rotating frame. However, this transformation Oz{^{t)) 
at the beginning and at the end of the evolution is close to the iden- 
tity operator, and taking it into account adds only a boundary effect, 
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which does not grow with the size of the time interval and is therefore 
neghgible. 

3.4 Berry phase under classical noise 

To find a dissipation-induced correction to the Berry phase we may 
use the same approach as in Section 2.3: first, we find the Lamb shift 
for a stationary field B and then evaluate the Berry phase using the 
relation (16). In this way we find the same expression (17) for the Berry 
phase. 

4. Conclusions 

In this paper we have derived expressions for the environment-induced 
correction to the Berry phase, for a spin coupled to an environment. On 
one hand, we presented a simple quantum-mechanical derivation for the 
case when the environment is treated as a separate quantum system. On 
the other hand, wc analyzed the case of a spin subject to a random clas- 
sical field. The quantum-mechanical derivation provides a result which 
is insensitive to the antisymmetric part of the random-field correlations. 
In other words, the results for the Lamb shift and the Berry phase are 
insensitive to whether the different-time values of the random-field oper- 
ator commute with each other or not. This observation gives rise to the 
expectation that for a random classical field, with the same noise power, 
one should obtain the same result. For the quantities at hand, our anal- 
ysis outlined above involving classical randomly fluctuating fields has 
confirmed this expectation. 

Furthermore, wc provided simple arguments, which allow one to un- 
derstand the contribution of fluctuations in various frequency ranges 
(below and above the Larmor frequency). 
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